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Abstract. Numerical results for the local field distributions of a family of Ising spin- 
glass models arc presented. In particular, the Edwards- Anderson model in dimensions 
two, three, and four is considered, as well as spin glasses with long-range power- 
law-modulated interactions that interpolate between a nearest-neighbour Edwards- 
Anderson system in one dimension and the infinite-range Sherrington-Kirkpatrick 
model. Remarkably, the local field distributions only depend weakly on the range 
of the interactions and the dimensionality, and show strong similarities except for near 
zero local field. 
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1. Introduction 

There has been interest in the distribution P{h, T) of local fields h at temperature 
T in spin glasses since the earliest days of their theoretical study [H [2]. Particularly 
infiuential was Thouless, Anderson and Palmer's [3] (TAP) self-consistent solution of 
P{h,T = 0) for the Sherrington-Kirkpatrick [1] (SK) infinite-ranged spin-glass model 
for which a mean-field theory is believed to be exact, albeit unusual and very subtle 
[SIE]. Since then there have been several further studies of P{h) for the SK model (see, 
for example, Refs. [7], [8], [9], and [10]) and the nature of the local field distributions is 
well understood. On the other hand, there has been little work on the study of P{h) 
for other spin-glass models. Most notably, few studies exist of P{h) for the finite-range 
canonical Edwards- Anderson [Til [12] (EA) Ising spin-glass model which is generally not 
exactly solvable. It remains controversial whether some of the specific subtleties of the 
SK model are applicable to the EA model and other more realistic models, and their 
relationship is far from clear. 

This paper studies and compares numerically the local field distributions at T = 
of Ising spin glasses with varying range-behaviour and spatial dimensionality, which 
are largely inaccessible to exact solution. We focus our discussion to Gaussian bond 
distributions of zero mean. A remarkable and somewhat surprising similarity is found 
across systems for which other aspects of the statistical physics state structure are 
believed to be different, but with systematic small differences near h = as the system 
interpolates between the limits of one-dimensional nearest neighbour and infinite-ranged 
mean field, at both extremes of which P{h = 0, T = 0) is zero in the thermodynamic 
limit. 

The paper is structured as follows: In Sec. [2] we introduce the models studied 
followed by brief descriptions of the numerical methods in Sec. [31 Our results are 
presented in Sec. jH followed by concluding remarks. 

2. Models 

Both the SK model and the EA models are characterized by the Hamiltonian 



where 5*^ G {±1}- The interactions Jij are chosen randomly and independently from 
a Gaussian distribution of zero mean and then quenched. In the SK model the sum 
over i and j extends over all sites and the variance of the distribution V{Jij) scales 
inversely with the system size as J'^ /N . Here we set J = 1 so that the spin glass 
onset transition is at = 1. For the EA model the sum over the indices i and j is 
restricted to nearest-neighbour pairs and the variance of the bonds Jij is independent 
of the number of spins A^ on a d-dimensional hyper-cubic lattice of size A^ = L*^; in 
this case, unlike for the SK model, the lattice dimension is relevant both qualitatively 
and quantitatively. There is universal agreement that the SK model exhibits a phase 




(1) 
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Figure 1. Schematic pliase diagram of the KAS model in the d-a plane following 
Ref. [15| . In this work we focus on = 1 which corresponds to the white horizontal 
arrow. For (t < 1/2 we expect infinite-range (IR) behaviour reminiscent of the SK 
model and where the energy needs to be rescaled as a power of the system size to 
avoid divergences, whereas for 1/2 < a < 2/3 we have mean- field (MF) behaviour 
corresponding to an effective space dimension doff > 6 (the thickened line separates 
mean- field from non- mean- field behaviour). For 2/3<(T<lwe have a long-range 
(LR"*") spin glass with a finite ordering temperature Tc, whereas for 1 ^ cr ^ 2 we 
have a long-range spin glass with Tc = (LR"). For cr > 2 the model displays short- 
range (SR) behaviour with a zero transition temperature |161 117j . Empirically, for 
0.5 < cr < 1 the d = 1 KAS model can be identified as corresponding to an EA system 
with effective dimension dos ~ 2/{2a — 1) [T3j. Figure adapted from Ref. |17j . 



transition as the temperature is reduced |T3] to a phase with an ultrametric hierarchy of 
metastable states and an associated mathematical feature of rephca-symmetry breaking 
of the overlap order parameter. For the finite-range EA model there is a spin-glass 
transition above the lower critical dimension, d > di, believed [H] to be di = 5/2, but 
the existence of an ultrametric hierarchy in short-range systems is controversial, not 
proven and disbelieved by many practitioners. 

In order to effectively interpolate between the physics of the SK and EA models 
and to probe their similarities and differences we also study a "tunable" model first 
introduced by Kotliar, Anderson and Stein (KAS) [16] and recently studied in detail 
by Katzgraber and Young [iTl [181 [13 [SO]- The model, which has helped elucidate 
many properties of spin glasses, is a long-range Ising spin glass with random power-law 
interactions. The Hamiltonian of the model is given by Eq. ([1]) but now with the sites 
i and j on a (i-dimensional lattice with periodic boundary conditions and the exchange 
interactions given by 




Here, is the separation of spins i and j, the e^j are chosen randomly and independently 
from a Gaussian distribution of zero mean and standard deviation unity, and c{a) is a 
constant. The KAS model is believed to interpolate between mean-field-like behaviour 
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Figure 2. Left panel: P{h) for h > for SK systems of sizes iV = 15, . . . , 1023 at 
T — 0. (Here, a bimodal bond distribution, J = -iil/^/N, has been used.) Clearly, the 
curves are well converged for the larger system sizes. The inset shows an enlargement 
of P{h) for h near zero. The slopes are all « 0.3, as indicated by the dashed line. 
Right panel: P{h = 0) 0.153(3)A^~^/^, i.e., the data extrapolate to zero in the 
thermodynamic limit. The low noise in the data suggest that errors are smaller than 
the symbols. 



for small a < aci{d), an intermediate long-range regime [aci{d) < a < crc2('^)]5 and 
a short-range regime [(Jc2{d) < cr < acsid) = oo]. Each of the latter two regimes are 
subdivided into ordering and non-ordering regimes depending upon the space dimension 
(higher dimensions favoring ordering); see Fig. [1] In the present work we shall consider 
only d = 1 for the KAS model, for which case the intermediate long-range region has 
a finite cooperative ordering temperature > for a < Uc but no finite-temperature 
ordering for a > a^, see Fig. [H To enforce periodic boundary conditions we place the 
spins on a ring (see Ref. [17] for details) and choose the geometric distance between the 
spins, i.e., = (N/tt) sin(7r|i — j\/N). We normalize the interactions to have = 1 
for all cr, i.e., 

c(a)-^ = $^r-2^ (3) 

For the short-range EA model we study hyper-cubic lattices in different space 
dimensions d. Again, in order to assist quantitative comparisons, in all cases we choose 
the exchange-scale normalization to yield the same mean-field transition temperature 

= 1. Hence for EA models we choose the variance of the exchange to be 1/z where 
z is the coordination number; for the hyper-cubic lattices that we study z = 2d. In the 
thermodynamic limit the SK model is equivalent to the infinite-dimensional EA model; 
hence its normalization with the variance of the exchange interactions scaling as 1/N. 
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3. Numerical procedures 

For the EA model we apply the Extremal Optimization (EO) heuristic as described in 
Ref. [21]. EO provides approximate ground states of spin glasses with high accuracy 
typically within 0{N^) update steps, at least for system sizes up to < 256 as studied 
here for the EA model. We read off and average the local fields for the presumed zero- 
temperature configuration found for each instance. For each reported system, between 
10^ — 10^ instances have been optimized, depending on system size. Since EO finds 
near-optima using a far-from-equilibrium dynamics, the explored configurations may 
possess a systematic bias. To check that this is not the case we have applied EO to 
reproduce P{h) for the SK model (with Jij e {±N-^/^}) for A^ < 1023, as studied 
in Ref. [22], finding no such biases. This is apparent from Fig. [21 in which familiar 
properties (discussed below) of P{h) in the SK model are reproduced, such as the linear 
slope for and the finite-size scaling of P(/i = 0) ~ 1/VN for N ^ oo. 

For the KAS model we use exchange (parallel tempering) Monte Carlo for the 
simulations [23]. To measure the local field distributions we compute 5000 disorder 
realizations for each system size A^ and value of a. The lowest temperature simulated 
is T = 0.05 which is close enough to T = such that for the system sizes studied we 
effectively probe the ground state of the system [2l] . We also note that earlier studies of 
the SK model [8] have shown that P{h = 0, T) = XT + 0{T^) with A ^ {27re)-^/^ ^ 0.25, 
giving a deviation of ~ 0.01 for the SK model at T = 0.05 compared with T = 0, which 
is indeed negligible compared to the values we find for small h for finite-range models. 

In this study we consider a = 0.00 (SK), 0.55 (MF), 0.75 (LR+), 0.83 (LR+), 1.00 
as well as 1.50 (LR°) and 2.00 (SR); see Fig. [1] for details. For all system sizes A^ 
and exponents a we study A^ = 29 replicas. For a < 1.0 we equilibrate the system for 
A'sw = 2^^ Monte Carlo sweeps for A^ < 128 and for 2^° Monte Carlo sweeps for A^ = 256. 
For 1.50 < cr < 2.0 we again take 2^^ Monte Carlo sweeps to equilibrate for A^ < 64 
but increase to 2^9 Monte Carlo sweeps for A^ = 128. Then in each case we measure 
the local field distributions for the same amount of Monte Carlo time. Equilibration is 
tested by comparing the energy calculated from the link overlap to the internal energy 
of the system calculated directly. Once both agree, the system is considered to be in 
thermal equilibrium. For details see Refs. [20] and |25j . 

4. Local field distributions 

The local field distribution is defined by 



where (■ ■ ■) indicates a thermodynamic average and [■ ■ -Jav denotes an average over the 
quenched disorder. 

The simplest nontrivial mean-field solution for P{h) for the SK or EA models, the 
replica-symmetric [HI HJ [2S] effective field approximation, yields P{h) self-consistently 
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through 

P{h) = ^=6-'''/^'^, q= [ dhP{h)tanh\ph), (5) 
V27rg J 

where f3 = 1/T is the inverse temperature. However, this approximation incorrectly 
yields a hole in P{h) as T ^ [26], i.e., P{h) = for all | h \< (2/7r)i/2. 

Indeed, Thouless, Anderson and Palmer [3] argued already in 1977 that for small 
fields h in the SK model P{h,T = 0) ~ 0.3|/i|, linearly in \h\ [27]. Later studies 
of the SK model, employing full replica symmetry breaking (FRSB) [2H], have borne 
out this linear form and the value of its slope (see, for example. Fig. [2] and Refs. [8], 
[To] , and [29]). Through accurate studies of a large sequence of finite-replica-symmetry 
breakings Refs. [9] and [10] have demonstrated that within replica theory the correct 
linear behaviour of P{h) for the SK model requires the full limit of infinite replica- 
symmetry breaking order; any finite-order truncation yielding a fictitious (if decreasing 
with RSB-order) hole in P{h) near h = 0. 

Another model that is exactly solvable is the one-dimensional nearest neighbour 
random-exchange Ising chain (limit of the KAS model for a ^ 1), the one- dimensional 
EA model [301 18]. It does not, however, have either frustration or a finite-temperature 
phase transition and, thus, no spin glass phase. P{h) is given by 



P{h) = J dJV{J)V{h + J) {1 - tanh[/3J] tanh [(3{h + J)]} , (6) 



with 



nJ) = ^e-^^ (7) 

71 



such that = 1/2. For P oo, this evaluates to 

P{h) = -\h\ rdxe-^'(2^'-2^+^), (8) 



" JO 

again giving P{h = 0) = and a linear small-/?, behaviour, analogously to the SK model, 
but with a slope almost twice as large, i.e., P{h) ~ 2\h\/7r. 

In higher dimensions the EA model is not exactly solvable. There is no finite- 
temperature spin-glass phase in d = 2, although there is such a phase [33] in d = 3 and 
greater. Yet, it is not clear that the glassy phase for d > 3 exhibits the characteristics 
of replica-symmetry breaking (e.g., ultrametricity) found in the SK model; at least up 
to an upper critical dimension believed to be du_cd = 6. We have performed a numerical 
simulation of P{h) for the EA model using the EO heuristic [21]. The results are 
exhibited in Figs. [31 [Since P{h) is symmetric, we only show plots for h > 0, i. e., 
P{h)dh = 1/2.] At first sight P{h) at T = shows very little variation between the 
dimensions d; the overall shape is quite similar to that for the SK model, see Fig. HI But, 
in fact, at closer detail there is a notable distinction for h ^ 0, where the behaviour 
for the EA model differs significantly from the SK model. Unlike for the SK model (see 
Fig. [2]), P{h = 0) appears to be finite in the thermodynamic limit of the EA model, 
as the plot near h = in the bottom panels of Figs. [3] indicate. We have extrapolated 
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Figure 3. P{h) for /i > for the EA spin glass on hyper-cubic lattices for d ~ 2, 3, 
and 4 for various side-lengths L for a Gaussian bond distribution with (J^) = l/{2d). 
Already beyond some small L, there is little distinction between the P{h) for different 
sizes, indicating only small corrections to scaling. This becomes even more apparent 
in the insets, showing an enlargement of the data near h = 0. P{h) in each d seems 
to converge to a finite value at /i = of P{h = 0) « 0.065 {d = 2), P{h = 0) « 0.055 
(d = 3), and P{h = 0) « 0.045 (d = 4), see Fig.O P{h) for small \h\ rises with a slope 
of a « 0.25 {d = 2), a sa 0.23 {d = 3), and a « 0.23 {d = 4), where P{h) - P{0) + a\h\. 
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Figure 4. P{h) for h > for the EA spin glass in d = 2, 3, and 4 dimensions and for 
the SK model, all at comparable system size N = 216 . . . 256 for bond distributions 
with (J^) = 1/z, where z is the connectivity of each spin: z = 2d for the EA model 
and z = — 1 for the SK model. Also plotted is the exact d = 1 result from Eq. ([5]). 
The plot highlights the overall similarity in P{h) for all models and space dimensions. 
Except for \h\ close to zero, the numerical data for the EA model in d = 2 - 4 overall 
seem to interpolate monotonically between the d = 1 result and the SK (d = oo) 
model. Yet, while P(0) = for both d = 1 and SK, for d = 2 - 4, P(0) is positive and 
monotonically decreasing, see Fig. [S] 
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Figure 5. Extrapolation of the values for P{Q) obtained for d = 2, 3, and 4 in Figs. [3] 
for various system sizes L. The extrapolation to the thermodynamic limit L oo 
proceeds on a scale of 1/L'^~y, where y for each d is the stiffness exponents discussed 
in Ref. [H]. As P{h) has units of inverse energy, these are the appropriate corrections 
to scaling for the EA model [311 [32]. Clearly, P{0) > for each d a.t L ^ oo, unlike 
for the SK model. The noise in the data suggests errors to be about double the size of 
the symbols. 
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h 



Figure 6. P{h) ior h > and different system sizes N for the KAS model in d = 1 
for different powers of the exponent a which change the effective space dimension. 
The data are for the lowest temperature simulated, T = 0.05. The insets show an 
enlargement of the area around /i = 0. P(/i = 0) clearly converges to a finite value 
for all a > 0.5 and shows an approximately linear behaviour [P{h) — P{0)] ~ 0.3|/i| 
for small \h\. Note that the corrections to scaling decrease considerably for larger a 
values. 
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Figure 7. Direct comparison of the local field distribution P{h) (at T = 0.05) as a 
function of the local field h for the KAS model and = 128 spins and different values 
of the exponent cr covering all possible universality classes from infinite range to short 
range with zero transition temperature. Also plotted are the exact results for d ~ 1 
(i.e., a ^ oo) from Eq. §^ and the SK data (i.e., a = 0) from Fig. [H (both at T = 0). 
As in Fig.m the data interpolate smoothly between both extremes, except at P(0). 

the values of -P(O) to infinite system size L in Fig. [51 For eacli d > 2, wlietlier above 
or below the lower critical dimension, P(0) quickly settles to a positive value between 
approximately 0.04 and 0.07. This value seems to decrease slowly with increasing space 
dimension d, consistent with P{0) = in the SK [d = oo) limit; see Fig. [13] below. 
P{h) — P{0) rises linearly for small h with a weakly d-dependent slope of a ~ 0.23 to 
0.25. 

We have performed a similar set of simulations for the d = 1 KAS model for a range 
of different a values covering different behaviours ranging from SK-like to finite-range 
non-ordering (Tc = 0) [IT]. In Fig. [6] we show the local field distributions for T = 0.05 
for different system sizes N. Each panel is for a different exponent a in Eq. ([2]). The 
similarity of all data sets for different a is clearly visible. The insets show always the 
area around /i = in detail. To further illustrate the similarities between the data sets, 
in Fig. [7] we show data for P{h) for = 128 and different exponents a covering all 
possible universality classes. The data for all a agree relatively well, with the data for 
cr > 1.0 showing a more pronounced peak and larger gap. Again, for the different values 
of cr > 0.5 studied, [P{h) — P(0)] ~ a\h\ for \h\ small with a pa 0.3 interpolating between 
the SK and the short-range one-dimensional result |34j . 

We have also extrapolated the data for P{h = 0,T, N) to T = (fits to a quadratic 
function for T < 0.3 with fitting probabilities [35] larger than ~ 0.3). A typical example 
of the extrapolation for a = 0.83 is shown in Fig. [8] (the behaviour of P{h) for different 
temperatures is shown in Fig. [9]). Since the difference between T = 0.05 and T = is 
minimal and because estimating the error bars for the extrapolated data is difficult, in 
Fig. [in] we show data for T = 0.05. 
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Figure 8. P{h = 0) as a function of temperature for a = 0.83 and N = 256. The 
data are well fit by a linear behaviour in T (with very small quadratic corrections) with 
slope « 0.25 for T < 0.3. Data for T = 0.05 are very close to T = (arrows), which is 
why in Fig. [TO] we extrapolate to = oo for T — 0.05 and not T = 0. Furthermore, the 
estimate of the error bars in the temperature extrapolation is difhcult. Inset: P(h = 0) 
for a range of a compared with the SK result (27re)^^/^r. 



A closer look at Figs. [7] and [TO] poses an interesting question: In the KAS model, 
P(0) appears to be an increasing function of a for all values studied so far, rising from 
P(0) = at a = for SK to P(0) ^ 0.08 at a = 2.0. Yet, the exact resuh for the 
one-dimensional Ising chain, corresponding to cr = oo, again has P{0) = 0. Hence either 
the apparent extrapolation to the limit ^ oo limit, shown in Fig. [TOl is incorrect 
and P(0) = for all a after all, or there has to be a maximum in P{0) at some finite 
o'max beyond which P{0) again descends to zero. To decide this question, we have done 
a more extended study of P{0) also for a > 2.0 which is shown in Fig. [Ill The results 
clearly show a well-defined maximum near (Tmax ~ 1-8, which is likely to persist in the 
thermodynamic limit. The resolution of this question — surprising in its own right — 
strengthens our belief that finite-size effects in our presentation are well under control. 
In all, the KAS model essentially reproduces the results found for the SK and canonical 
EA models for cr < 1 where it corresponds to possible physical dimensions and for 
larger a [where the formula dcs ~ 2/ {2a — 1) [I3] becomes inappropriate] it continues 
smoothly towards the nearest-neighbour d = 1 EA limit. The fact that P{0) peaks for 
the KAS model at a value of a intermediate between that corresponding to des = 2 and 
des = 1 suggests that if one could continue the EA model off integer dimensions, there 
might exist a dimension dma.x, likely between d = 1 and d = 2 at which P(0) would 
peak. The slope of [P{h) — P(0)] near h = is however found to vary monotonically 
with a between the SK [a = 0) and nearest neighbour EA {a = oo) limits. 

The local field distribution P{h) for a disordered Ising spin system with Gaussian 
distributed random exchange disorder is thus seen to be remarkably robustly linear in 
\h\ at low h with a coefficient in the range a = 0.25 - 0.35, irrespectively of whether the 
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Figure 9. P{h) for cr = 0.83 and N = 128 for different temperatures ranging from 
T = 0.05 to temperatures well above the transition temperature T^. For low enough T 
the local-field distribution shows a dip around h = which is "filled in" for increasing 
temperatures. For T 3> Tc we obtain a Gaussian distribution. 




Figure 10. Extrapolation of the values for P(0) obtained for the KAS model as a 
function of 1/Vn for T — 0.05 and different values of the exponent a. The data for 
a = 0.00 decay with a power-law similar to the SK model. For larger values of a the 
data indicate a saturation to a finite value of P{0) for N ^ oo. Larger system sizes 
would be needed to clearly differentiate the different scenarios beyond a qualitative 
comparison. The noise in the data suggests an error of about 5 times the symbol 
size. For comparison we also show the data for SK at T = from Fig. [21 note that 
the shift compared with the cr = KAS results is due to the finite-temperature shift 
(2^e)-i/2r. 



system exhibits spin-glass order and independently of whether such spin-glass order is 
replica-symmetry- broken mean field or not. 

It is also interesting to compare with studies of a non-equilibrated SK model 
[36l [371 [38] . In these studies, spins are initially randomized and then exposed to purely 
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Figure 11. Plot of P{0) as a function of cr in the KAS model. For each data point in 
= 32 and 64 averaged P{h) over at least 10^ instances, 2 10'' instances at iV = 128, 
and at least 10'^ for N = 256. (The data for N = 256 clearly show a systematic bias 
only near the maximum, indicative that those instances are hardest to optimize. Note 
the complete lack of residual finite size effects for a > 2.0). To obtain smooth estimates 
of P(0), independent of the width of histograms chosen to measure P{h), we made a 
linear fit to the respective P{h) for h < 0.3 to obtain each data point. Inset: Plot of 
the slope a = dP{h)/dh\h-,a of P{h) near /i = as a function of a, extracted from 
a fit to the same data. Only data for N = 256 are shown, as there are hardly any 
finite size effects. 2/7r (solid line) represents the slope for the one-dimensional result 
in Eq. ([5]) obtained for cr ^ oo. 



relaxational single-spin dynamics, i.e., dynamics in which spins are chosen randomly 
and flipped if and only if such a flip would reduce the energy, until a metastable state 
is reached. The distribution P{h) averaged over only those metastable states reached 
by this procedure is again linear in \h\ for small h with a slope that again appears to 
be of order a ~ 0.3. The interest of this observation in the present context is that the 
metastable states reached in this dynamical procedure are not the ground state. Indeed 
they are significantly above the ground state, with the average relaxational energy 
variously reported as E^c\a.y./N ~ —0.7 [37], -0.715 [36], and —0.73 [38j, much higher 
than the ground state energy per spin Egg/N = —0.7632 [TD], while the corresponding 
P{h, -Erciax) averaged over all metastable states has a finite value at h = [361 137] . 

We have conducted a comparable study of P{h) over metastable states reached by 
rapid-quench for the EA model in d = 2, 3, 4, and 7, and SK, as shown in Fig. [121 As for 
the results for the SK model, there are significant (qualitative) similarities between the 
T = equilibrium results oiP{h) in Figs.[3land those obtained by simple relaxation. But 
especially the behaviour near h = deviates quantitatively, with P{0) distinctly larger 
and the initial slope significantly smaller. Yet, both seem to tend smoothly towards 
the corresponding SK result P{0) = (see Fig. Hand Refs. [36] and [37], respectively) 
for d — * oo, as is demonstrated in Fig. [131 Furthermore, for small \h\ the local field 
distribution shows a linear behaviour as found in Ref. [39] via an intrinsically far-from- 
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Figure 12. Local field distribution P{h) for the EA model in c? = 2, 3, and 4, obtained 
after purely rclaxational single-spin dynamics on 10* instances for each L and d. The 
same system sizes L as in Figs. [3] have been used here but finite-size effects are well 
below the statistical noise. We have added results ior d ~ 7 {L ~ 3) to show that 
there are no drastic changes above the upper critical dimension ducd — 6. The dashed 
line merely serves to guide the eye. Inset: P{h) for relaxational single-spin dynamics 
in SK. 
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Figure 13. Extrapolation of P{Q) towards the SK-limit, d oo. Shown arc the 
values of P(0) obtained from the L oo extrapolation in Fig. [5] for ground states 
of the EA model (full circles) and those obtained by a relaxation process of the data 
shown in Fig. 1121 ffuU diamonds), which are approximately L-indepcndent. Added are 
the data points for SK for ground states from Fig. [2] (open circles) and for relaxation 
from Fig. [12] (open diamonds). As the EA spin glass approaches the SK model by 
setting 2(i = A^ — 1 — *■ oo, the appropriate scale here is (2(i)~^/^ ~ N~^^'^, according to 
the right panel of Fig. O Asymptotic slopes for the SK data are 0.153(3) for ground 
states and approximately 0.61 for relaxation (see also Ref. [3H])- For both sets of data, 
it seems that limd^oo^'(O) 0, consistent with the respective SK results. Similar to 
the case of tr > 0.5 in Fig. [TTl there is no indication here that P(0) = for any d < oo. 
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equilibrium simulation along the hysteresis loop of the model for finite external fields. 
Finally, as in Refs. [36] and [37], the metastable states obtained under relaxation in 
the EA model are substantially more energetic than the ground states. For instance, 
in c? = 3 relaxation gives a normal distribution of states of mean ^'reiax/-^ ~ —1.4 
and a deviation of approximately 0.1, whereas the ground states found at L = 6 are 
centered just above Egs/N ^ —1.7, with a much narrower deviation and close to the 
thermodjTiamic value of —1.700(1). [IQ] 

Finally in this section, let us return to Fig. [8] now concentrating on the gross 
T-dependence. This figure demonstrates another robust feature, the linearity of 
P{h = 0,T) with T up to the mean-field transition temperature Tmf = 1; this result 
previously observed in the SK model [8] is seen from the main figure to hold well also 
for the KAS model with a = 0.83, even though the actual transition temperature in this 
case is much smaller than the mean-field temperature (indeed it is closer to Tc = 0.45 
|20j). As shown in the inset, a similar behaviour is found for other values of a (up to 
the order of the peak in Fig. [TTl) . the higher of which have no finite-temperature phase 
transition. Thus again the general shape of P{h = 0,T) — P{h = 0,T = 0) is rather 
robust against whether the system is one which orders or not, has RSB or not. 

5. Conclusions 

We have presented results for the local field distributions of different spin-glass models 
in different space dimensions. These include both cases where there is generally believed 
to be a finite-temperature phase transition and others where no finite temperature 
transition occurs. They also include cases where the subtleties of replica symmetry 
breaking are believed to operate and others where they do not occur or are in question. 
Our results show that the distributions of local field are qualitatively very similar for 
all models with Gaussian-distributed interactions; the data for P{\h\,T = 0) decay 
exponentially for large fields \h\ and show a linear behaviour for \h\ close to zero. The 
distributions only show differences in behaviour near T = and h = where for the 
infinite-range SK model P{h = 0,T = 0,iV) ~ iV-^/^ ^ g as ^ cx), while all other 
finite-dimensional Edwards- Anderson spin-glass models studied with space dimensions 
d > 1 P{h = 0, A^) seems to tend to a constant in the thermodynamic limit. For finite d 
the thermodynamic value of -P(O) scales as d~^^'^ for low d with a coefficient very close 
to that for the SK model, taking 2d = N — 1 to match coordination number for spins 
in EA and SK, see Fig. [131 This suggests that (i~^/^-scaling applies for all d |41] . 

These observations are mirrored by simulations of a one-dimensional Ising model 
with random power-law interactions: for the regime of the power-law exponent a which 
correspond to an infinite-ranged system, P{h = 0) decays with an inverse power of the 
system size, whereas for all other universality classes the distributions tend to a constant 
in the thermodynamic limit, rising from zero in the limits of both SK a = and the 
unfrustrated nearest-neighbour a = oo, with a maximum for a near 2. Furthermore, for 
all a the P{h = 0,T) are all close to the same linear-T behaviour for T < Tmf = 1- 
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Qualitatively similar, but quantitatively different, small-/i P{h) behaviour is also 
found for systems that are quenched from random starts (to states that are not 
thermodynamically equilibrated). 

Our study has been concerned with the case of Gaussian exchange distribution. 
For problems with discrete distributions, such as for Jij randomly ±J one expects 
P{h = 0,T = 0) to be nonzero for finite-range systems [8]. It is, however, surprising 
how small the observed P{h = 0,T = 0) are for Gaussian exchange-disorder without 
being zero. This effectively counsels against associating small -P(O) and simple small-h 
slopes with the subtleties of RSB or finite-temperature spin glass transitions [i2] . 

Finally, we should caution that the simulations were performed for finite system 
sizes and although corrections to scaling seem to be very small, a change in behaviour 
at larger system sizes cannot be ruled out completely. 
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